Route to supersolidity for the extended Bose-Hubbard model 
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We use the Gutzwiller ansatz and analyze the phase diagram of the extended Bose-Hubbard Hamiltonian with 
on-site {U) and nearest-neighbor iV) repulsions. For d-dimensional hypercubic lattices, when 2dV < U, it is 
well-known that the ground state alternates between the charge-density-wave (CDW) and Mott insulators, and 
the supersolid (SS) phase occupies small regions around the CDW insulators. However, when 2dV > U, in this 
paper, we show that the ground state has only CDW insulators, and more importantly, the SS phase occupies 
a much larger region in the phase diagram, existing up to very large hopping values which could be orders of 
magnitude higher than that of the well-known case. In particular, the SS-superfluid phase boundary increases 
linearly as a function of hopping when 2dV > 1.5U, for which the prospects of observing the SS phase with 
dipolar Bose gases loaded into optical lattices is much higher. 

PACS numbers: 03.75.-b, 67.80.kb, 67.85.-d, 67.85.Hj 



Can a solid be superfluid? The so-called supersolid phase is 
characterized by the simultaneous existence of diagonal (crys- 
talline) and off-diagonal (superfluid) long-range orders. Al- 
though this intriguing possibility was suggested a long time 
ago in the context of solid "^He UJ, and in spite of numerous 
attempts over the past decades, a convincingexperimental ev- 
idence for its existence is yet to be found [2]. On one hand, 
there is still some controversy in the condensed matter litera- 
ture around the recent reports that the theoretically predicted 
nonclassical rotational inertia was found in solid ^He with the 
torsional oscillator experiments [3]. On the other hand, there 
is strong theoretical evidence that the situation in lattice mod- 
els is promising fllH, which could be advanced with ultra- 
cold quantum gases loaded into optical lattices 

Possibly the simplest models that show SS behavior are the 
extended-type Bose-Hubbard ones with on-site and nearest- 
neighbor (NN) repulsions. Since these models can be nat- 
urally realized with dipolar bosons i-e. bosonic atoms 
or molecules with permanent or induced magnetic or electric 
dipole moments, the ground-state phase diagrams of various 
extended models have already been studied. For instance, 
the existence and stability of SS phases have been demon- 
strated via the Gutzwiller ansatz S 01 and decoupling 
mean-field IIC app roaches, and numerically exact quantum 
Monte Carlo js. flltl techniques. The SS is known to be a very 
fragile phase, and one of the major obstacles in creating and 
observing it is its very small existence region and low critical 
temperatures. In a very recent proposal 1 12], it has been sug- 
gested that one way of overcoming these difficulties is to load 
high occupancies of bosons into optical lattices. Such systems 
are well-described by the quantum rotor model, for which the 
mean-field calculation gives linear dependence between the 
critical hopping (tc) and occupancy. 

In this paper, we use the Gutzwiller ansatz and mean-field 
decoupling to analyze the phase diagram of the extended 
Bose-Hubbard Hamiltonian with on-site (U) and isotropic NN 
(V) repulsions. For c?-dimensional hypercubic lattices, be- 
yond the critical threshold 2dV > U, we show that the SS 
phase occupies a much larger region in the phase diagram, ex- 



isting up to very large hopping values of the order 2dtc > U. 
In particular, the SS-superfluid phase boundary increases lin- 
early as a function of hopping when 2dV > 1.5U. This 
must me contrasted with the below threshold {2dV < U) 
case, for which it is well-known that the SS phase occupies 
small regions around the CDW insulators, existing only up to 
2dtc < QAU . Therefore, we show that the prospects of ob- 
serving the SS phase is much higher when 2dV > U. We also 
argue that our results for two dimensions is directly applicable 
to the quasi-two-dimensional dipolar Bose gases, for which 
the condition AV > U could be easily achieved by tuning the 
s-wave scattering length using the Feshbach resonances. 

Hamiltonian: To obtain these results, we use the extended 
Bose-Hubbard Hamiltonian with an isotropic NN repulsion 



H = - tJ2iblbj + blh) + - 1) 



(1) 



where t is the tunneling (or hopping) amplitude between NN 
sites i and j, b\ (bi) is the boson creation (annihilation) oper- 
ator at site i, rii = b\bi is the boson number operator, and /i 
is the chemical potential. As it turns out, the phase diagram 
of this Hamiltonian depends strongly on the relative strength 
U — zV, where z = 2o? is the lattice coordination number, 
leading to two important cases. 

(I) Weak-NN-coupling (zV <U): When V ^ 0, the ground 
state has two types of insulating phases H, 15]. The first one 
is the Mott insulator where, similar to the usual Bose-Hubbard 
model, the average boson occupancy is the same for every lat- 
tice site, i.e. (n^) = rtQ. Here, (...) is the thermal average, 
and no is chosen to minimize the ground-state energy for a 
given fi. The second one is the CDW insulator which has 
crystalline order in the form of staggered average occupan- 
cies. To describe the CDW insulator, it is convenient to split 
the entire lattice into two sublattices (e.g. A and B) such that 
the NN sites belong to a different sublattice, i.e. (n^) = ua 
and (rij) = ub for We assume the occupancies are 
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such that UA > nB, and the case with ua = ns corresponds 
to the Mott insulator 

In the atomic {t = 0) limit, it turns out that the chemical 
potential width of all CDW and Mott insulators are zV and 
U, respectively, and that the ground state alternates between 
the CDW and Mott phases as a function of increasing /i. For 
instance, the ground state is a {ua — 0, ub = 0) vacuum for 
/i < 0; a (1, 0) CDW insulator for < /i < zV; a (1, 1) Mott 
insulator foTzV<fi<U + zV; a (2, 1) CDW insulator for 
U + zV<ii<U + 2zV- a (2, 2) Mott insulator for U + 
2zV < fi < 2U + 2zV, and so on. As t increases, the range 
of n about which the ground state is insulating decreases, and 
the CDW and Mott insulators disappear at a critical value of 
t, beyond which the system becomes compressible (SF or SS) 
as shown in Fig.|2a). 

(II) Strong-NN-coupling (zV > U): In contrast to the 
well-known weak-NN-coupling, the strong-NN-coupling of 
this model has not been studied much in the literature, which 
is the main topic of this paper. When zV equals exactly to 
U, it is easy to check at least in the atomic limit that the 
(no + l,no) CDW insulator becomes degenerate in energy 
with the (2no + 1,0) CDW insulator, and the (no, no) Mott 
insulator becomes degenerate with the (2no, 0) CDW insula- 
tor. This indicates that, beyond the critical zV = U threshold, 
both the CDW and Mott insulators that are found in the weak- 
NN-coupling are unstable against formation of new CDW in- 
sulators, leading to a phase diagram which has a very different 
qualitative structure. 

In this paper, we show that the ground state has only CDW- 
type insulating phases in the strong-NN-coupling, the chem- 
ical potential width of all are U in the atomic limit. For in- 
stance, the ground state is a (0, 0) vacuum for /i < 0, a (1, 0) 
CDW insulator for < ^ < C/, a (2, 0) CDW insulator for 
U < n < 2U, a (3,0) CDW insulator for 2U < < W, 
and so on. As t increases, the range of /i about which the 
ground state is insulating decreases, and the CDW insulators 
disappear at a critical hopping, beyond which the system be- 
comes a SS, as shown in Figs.|2b-f). Most importantly, unlike 
the weak-NN-coupling where the SS phase occupies small re- 
gions in the phase diagram around the CDW insulators, we 
show that the SS phase occupies a much larger region in this 
case, existing up to very large hopping values which could 
be orders of magnitude higher than that of the weak-NN- 
coupling. To obtain these results, we solve the Schrodinger 
equation for the Gutzwiller ansatz, as discussed next. 

Gutzwiller ansatv This ansatz has been frequently used to 
approximate the many -body wave functions of Bose-Hubbard 
Hamiltonians 141^^1. It can be written as 



choose m„ 



50. The normalization of the wave function 



n E/^^ 



(2) 



where \i,m) represents the Fock state of m bosons occupying 
the site i, and fi^m is the probability amplitude of its occupa- 
tion. Here, m = 0, 1, . . . , ra^^j^, where m^ax is the maximum 
number of bosons allowed in the numerics which we typically 



(■0I-0) requires = 1 for each i. 

Within this ansatz, the superfluid order parameter 
{il)\bi\il]) is determined by 



= ^ l/*„J^,„+l. 



(3) 



This complex parameter describes the state of the system at 
site i: while it vanishes for the CDW and Mott insulators, it 
is finite for the SF and SS ground states. Therefore, 1/)^ -> 
0+ signals the phase boundary between an insulating and a 
compressible phase. Similarly, the average occupancy n^ — 
{4'\h\bi\%p) is determined by 



(4) 



In Eqs. (O and (|4|i, the probability amplitudes are obtained by 
solving the Schrodinger equation, (7/;|i/|V-') = ih{ilj\d\'4)) / dr, 
with = /o,^,me-"°■■"/^ This leads to 



-t (0iV^/o,'i,m-l + (?!)*Vni+ l/o,i,m+l) 



—m{m — 1) + Vmrii — /im 



/o 



(5) 



nj sum over sites j 



where (j}i = 2l(j). 4>] and fn 
neighboring to site i. 

For the ground state, first we need the minimal eigenvalue 
eo,j and the elements fo,i.m of the corresponding eigenvec- 
tor, and then use them in Eqs. (|3]l and (|4]l to solve for (pi and 
Tii self-consistently. Note that, in uniform systems with two 
sublattices discussed in this paper, the state of the whole sys- 
tem is sufficiently described by two order parameters: (pA and 
(j)B for the sublattices A and B, respectively. When this is the 
case, note also that eo,i = eo,A, fo,i,m = /o,A,m, (t>i = z(I)b, 
Ui — ua and n, = zub for i E A sublattice, and eo,i = eo,B, 

/o,i,m = fo,B,m, (pi = z4>A, Ui = UB and Hi = ZHA foV 

i e B sublattice. To support our Gutzwiller ansatz calcula- 
tions, next we examine the mean-field theory, which provides 
an analytical expression for the phase boundary between the 
insulating and compressible phases. 

Mean-field decoupling approximation: In constructing the 
mean-field theory, one first defines the SF order parameter 
4'i — {bi), and then replaces the operator bi with (f>i + 5bi 
in the hopping terms of Eq. ([T]l. This approximation decou- 
ples the two-particle hopping terms into single-particle ones, 
and the resultant mean-field Hamiltonian can be solved via 
exact diagonalization in a power series of (pi. Performing a 
second-order perturbation theory in (pi around the insulators, 
and following the usual Landau procedure for second-order 
phase transitions, i.e. minimizing the energy as a function of 
(pi, we eventually arrive at the condition 



(pi = (l>it 



1 



Urii + Vrii — p U{ni — 1) + Vrii — p 



(6) 
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where the definitions of <f>i and fii are the same as in Eq. (|5]). 
It is known that the results of the mean-field theory coincide 
with those of the Gutzwiller ansatz, and that they both become 
exact when d ^ 1 11311 . We emphasize that the mean-field 
calculations gives a good qualitative description of the system, 
and it becomes progressively accurate as the dimensionality 
and/or the occupancy increases. 

For uniform systems with two sublattices, Eq. ^ gives cou- 
pled equations for (j>A and (fiB, which can be solved to obtain 
the phase boundary between the insulating (CDW or Mott) 
and compressible (SF or SS) phases. Since (pAjfpB ^ 0+ 
near these boundaries, Eq. (|6]l can be satisfied only if Ij^ 
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UnB + zVnA — ^-t U{nB — I) + zVuA — fJ- 

(7) 

which gives a quartic equation for fi. An alternative way of de- 
riving this equation is the random-phase approximation lIlSll . 
Since a simple closed form analytic solution for fi is not pos- 
sible, we solve Eq. (|7]) for each of the insulating lobes sepa- 
rately. Having discussed the details of the Gutzwiller ansatz 
and mean-field approximation, we are ready to discuss the 
phase diagrams. 

Phase diagrams: We solve Eqs. ([3]|5]l self-consistently for 
the order parameters (i.e. ipA and (ps) and average occupan- 
cies (i.e. UA and ub), and use them to construct the phase 
diagram of the system. The CDW and Mott insulators are 
characterized by cj)A ^ (pB ~ 0, and ua ^ hb and ua — tlb, 
respectively. However, the SF and SS phases are characterized 
by (/>A = '/'B 7^ and cpA ^ <Pb, respectively, which naturally 
leads to ua = tlb in the SF and ua ^ tlb in the SS phase. 
In this paper, we choose cpA and (pB to be real, since we are 
only interested in the CDW-SS and MI-SF phase boundaries 
which are determined by {(pA, 0s} 7^ 0, and the SS-SF phase 
boundaries which are determined hy (pA = <pB 0- 

Therefore, it is sufficient to look at the relative order pa- 
rameter <pA — (pB and relative average occupancy n^i — to 
distinguish between various phases. For instance, in Fig.lT] we 
show the colored maps of (pA ~ <Pb and ua — tlb as a function 
of the chemical potential /i and hopping t for zV — 1.15U. 
The three dark lobes shown in Fig.lTfa), where (pA = <Pb = 0, 
correspond to (1,0), (2,0) and (3,0) CDW insulators (from 
bottom to top), the occupancies of which are clearly seen in 
Fig lfb). Here, the SF phase occupies the dark region that is 
common in both Figs. [Tfa) andfllb). It is also clear that the 
SS phase, where (pA — <pB and ua ~ riB ^ shown 
with bright colors in Figs. [Ha) anddfb), respectively, is sand- 
wiched between the CDW insulators from the left and SF 
phase from the right. Note that, in the SS phase, the super- 
fluid order parameter is larger on the sublattice with higher 
occupancy, and both the crystalline and superfluid orders are 
primarily on the same sublattice (13 • This is because the par- 
ticle and hole excitation energies are higher on the sublattice 
with lower occupancy when zV > U . 
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FIG. 1. (Color online) The colored maps of the (a) relative order 
parameter 0a — <^s and (b) relative average occupancy ua — ub are 
shown as a function of the chemical potential /i and hopping t when 
zV = 1.15C/. 



In Fig. |2] we repeat this analysis for a number of NN re- 
pulsions, and plot the phase diagrams as a function of /i 
and t. In these figures, the red continuous and blue dotted 
lines are obtained from the Gutzwiller ansatz calculations, and 
the black dashed lines are obtained from Eq. (|7]i. Note that 
both methods are in complete agreement, i.e. on top of each 
other, for the phase boundary between the insulating (CDW 
or Mott) and compressible (SF or SS) phases, which supports 
our Gutzwiller ansatz calculations. 

A typical weak-NN-coupling phase diagram is shown in 
Fig.|2a) for zV — 0.9C/. As discussed above, the ground state 
alternates between the CDW and Mott insulators as a function 
of increasing /i, and the SS phase occupies only small regions 
around the CDW insulators. Unlike the weak-NN-coupling 
phase diagrams. Figs. [2b-f) show that the ground state has 
only CDW-type insulating phases in the strong-NN-coupling. 
In particular, the chemical potential width of all CDW insula- 
tors is U in the atomic limit, and the ground state is a (?iOi 0) 
CDW insulator for (tiq — 1)U < n < tiqU . As t increases, 
the range of /i about which the ground state is insulating de- 
creases, and the CDW insulators disappear at a critical value 
of tc, beyond which the system becomes a SS. Note that, ex- 
cept for the weak-NN-coupling phase diagram, all of the tc 
values are comparable to each other for a fixed zV. 

The main result of this paper is shown in Figs. |2tc-f). It 
is clearly seen that the SS phase occupies a much larger re- 
gion when zV > I.IU, existing up to very large hopping 
values which could be orders of magnitude higher than that of 
the weak-NN-coupling. This result is intuitive given that the 
ground state has only {uq , 0) CDW insulators whose sizes are 
comparable to each other, and that the CDW modulations be- 
come stronger as no (or /i) increases. In fact, the SS-SF phase 
boundary becomes linear in /i and zt when zV > 1.5U. Our 
numerical calculations suggest that the slope of this line is ap- 
proximately given by dfj,/d{zt) = 2/{zV/U — 1). A similar 
linear dependence between the particle density and hopping 
has recently been found for the SS-SF phase transition bound- 
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FIG. 2. (Color online) The ground-state phase diagrams are shown 
as a function of the chemical potential and hopping t for (a) 
zV = 0.9(7, (b) zV = 1.1(7, (c) zV = 1.25(7, (d) zV = 1.517, (e) 
zV = 2(7, and (f) zl/ = 4(7. The CDW and Mott insulators are indi- 
cated with their sublattice occupancies (n^, ns). Here, the red con- 
tinuous and blue dotted lines are obtained from the Gutzwiller ansatz 
calculations, and the black dashed lines are obtained from Eq. Q. 



NN repulsion, where Cm — fJ-aP^ (or p'^ /eq) for particles with 
permanent magnetic (or electric) dipole moment p. The ratio 
U/{W) = 7r2^/2¥s3/4/(i2e^d), where em = add/a^ and 
add = mCdi/ {12Tih?) is the dipolar length scale, determines 
the critical threshold {AV > U). For instance, p — 6^b and 
add = 16ao for the ^^Cr atoms ifivll . where ^ib (ao) is the 
Bohr magneton (radius), and the condition is < Sao/s'^^^. 
However, p — 0.6 Debye and add = 2 x lO'^ao for the KRb 
molecules ITtIi . and the condition is < lO'^ao/s'^/*. Since 
tts ^ lOOan for Cr atoms, needs to be tuned via a Feshbach 
resonance jlSIl in order to achieve the critical threshold. 

In addition, we note that it is easy to extract the finite-size 
effects of an external trapping potential (e.g. present in atomic 
systems) from Fig. |2] For instance, within the local-density 
approximation, if the center of the trap is a (2, 0) CDW indu- 
lator (say zt < 0.6(7), the system is expected to go through 
first a SS, then a (1, 0) CDW indulator, before becoming a SF 
as a function of the radial distance towards the edge of the 
trap. On the other hand, if the center of the trap is a SS (say 
zt > 0.6(7), the system is expected to become a SF beyond a 
critical radius, without any intermediate phase. 

Conclusions: To summarize, beyond the critical threshold 
zV > U, we showed that the SS phase occupies a much larger 
region in the phase diagram, existing up to very large hopping 
values of the order zt > U, and that the SS-SF phase bound- 
ary increases linearly as a function of t when zV > 1.5(7. 
Therefore, our results suggest that the prospects of observing 
the SS phase is much higher when zV > (7, which could 
be easily achieved with quasi-two-dimensional dipolar Bose 
gases loaded into optical lattices, by tuning the s-wave scat- 
tering length via currently available Feshbach techniques flst]. 
We believe our mean-field treatment captures the qualitative 
physics right, and that this work will motivate further quantum 
Monte Carlo calculations in the strong-NN-coupling regime 
for more accurate phase diagrams. 

The author thanks J. K. Freericks for comments. This work 
is financially supported by the Marie Curie International Rein- 
tegration (Grant No. FP7-PEOPLEIRG-2010-268239) and 
the Scientific and Technological Research Council of Turkey 
(Career Grant No. TUBiTAK-3501-110T839). 



ary in the case of quantum rotor model 11211 . when the average 
occupation is much higher than unity. 

Experimental realization: Here, we argue that our results 
for two dimensions (z = 4) is directly applicable to the 
quasi-two-dimensional dipolar Bose gases. For the optical 
lattice potential VoL{r) = Vb[sin^(fcx) + sin^(fcj/)], where 
k = 27r/A is the wave vector and £ — A/2 the lattice 
spacing, the on-site interaction depends on the s-wave scat- 
tering length as via f^U = y^kasErS^^\ Typically, 
s — Vo/Er ~ 10, where Er — h?k^/{2m) is the recoil en- 
ergy and m is the particle mass. Assuming all of the dipoles 
are polarized along the z direction, the dipole-dipole inter- 
action becomes isotropic, leading to ^ = Cm/{4:Tt£^) for the 
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